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We present a study of the dijet invariant mass spectrum in events with two jets produced in 
association with a W boson in data corresponding to an integrated luminosity of 4.3 fb _1 collected 
with the DO detector at y/s = 1.96 TeV. We find no evidence for anomalous resonant dijet production 
and derive upper limits on the production cross section of an anomalous dijet resonance recently 
reported by the CDF Collaboration, investigating the range of dijet invariant mass from 110 to 
170 GeV/c 2 . The probability of the DO data being consistent with the presence of a dijet resonance 
with 4 pb production cross section at 145 GeV/c 2 is 8 x 10~ 6 . 
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The CDF Collaboration at the Fermilab Tevatron pp 
collider recently reported a study of the dijet invari- 
ant mass (Mjj) spectrum in associated production with 
W -> tv (£ = e or ,u) at s/s = 1.96 TeV with an in- 
tegrated luminosity of 4.3 fb _1 1]. In that paper they 
present evidence for an excess of events corresponding to 
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3.2 standard deviations (s.d.) above the background ex- 
pectation, centered at Mjj = 144 ± 5 GeV/c 2 >]. The 
CDF authors model this excess using a Gaussian peak 
with a width corresponding to an expected experimental 
Mjj resolution for the CDF detector @| of 14.3 GeV/c 2 
and further estimate the acceptance and selection effi- 
ciencies by simulating associated W + Higgs boson (H) 
production in the decay mode H — »• bb and with a mass 
M H = 150 GeV/c 2 . Assuming the excess is caused by 
a particle X with B(X -> jj) = 1, the CDF Collabo- 
ration reports an estimated production cross section of 
a(pp -> WX) 4 pb. 

Using 5.3 fb _1 of integrated luminosity, the DO Col- 
laboration has previously set limits on resonant bb pro- 
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duction in association with a W boson in dedicated 
searches for standard model (SM) Higgs bosons in the 
WH -> Ivbb channel @|. The DO Collaboration reported 
upper limits on a(pp — > lUiJ) x ,6(iJ — > 66) ranging 
from approximately 0.62 pb for Mh = 100 GeV/c 2 to 
0.33 pb for M H = 150 GeV/c 2 . The CDF Collaboration 
has performed a similar search using 2.7 fb _1 of inte- 
grated luminosity and reported no excess of events 
Furthermore, the DO Collaboration has not observed a 
significant excess of associated W boson and dijet pro- 
duction in analyses of either WW/WZ — > Ivjj [5| or 
H -> WW -> Ivjj @ using 1.1 ftr 1 and 5.4 fb" 1 of 
integrated luminosity, respectively. 

In this Letter we report a study of associated W(—> iv) 
and dijet production using data corresponding to 4.3 fb -1 
of integrated luminosity collected with the DO detector [7] 
at y/s = 1.96 TeV at the Fermilab Tevatron pp Collider. 
The CDF study of this production process uses the same 
integrated luminosity. We investigate the dijet invari- 
ant mass range from 110 to 170 GeV/c 2 for evidence of 
anomalous dijet production. 

To select W{— > Iv) + jj candidate events, we impose 
similar selection criteria to those used in the CDF analy- 
sis: a single reconstructed lepton (electron or muon) with 
transverse momentum pr > 20 GeV/c and pseudorapid- 
ity [|[ \r)\ < 1.0; missing transverse energy ]£t > 25 GeV; 
two jets reconstructed using a jet cone algorithm Q with 
a cone of radius A1Z = 0.5 that satisfy pr > 30 GeV/c 
and \r]\ < 2.5, while vetoing events with additional jets 
with pr > 30 GeV/c. The separation between the two 
jets must be | A^Qeti, jet2)| < 2.5, and the azimuthal 
separation between the most energetic jet and the di- 
rection of the TpT must satisfy A</>(jet, $t) > 0.4. The 
transverse momentum of the dijet system is required to 
be pr(jj) > 40 GeV/c. To reduce the background from 
processes that do not contain W — > Iv decays, we require 
a transverse mass of Mf? > 30 GeV/c 2 . In addi- 
tion, we restrict Mj. v < 200 GeV/c 2 to suppress muon 
candidates with poorly measured momenta. Candidate 
events in the electron channel are required to satisfy a 
single electron trigger or a trigger requiring electrons and 
jets, which results in a combined trigger efficiency for 
the evjj selection of (98l 3 )%. A suite of triggers in the 
muon channel achieves a trigger efficiency of (95±5)% for 
the \iv jj selection. Lepton candidates must be spatially 
matched to a track that originates from the pp interac- 
tion vertex and they must be isolated from other energy 
depositions in the calorimeter and other tracks in the 
central tracking detector. 

Most background processes are modeled using Monte 
Carlo (MC) simulation as in the CDF analysis. Di- 
boson contributions (WW, WZ, ZZ) are generated 
with pythia [ll| using CTEQ6L1 parton distribution 
functions (PDF) [T3|. The fixed-order matrix element 
(FOME) generator ALPGEN with CTEQ6L1 PDF 
is used to generate VF+jets, Z+jets, and it events. The 



FOME generator COMPHEP [14j is used to produce sin- 
gle top-quark MC samples with CTEQ6M PDF. Both 
ALPGEN and COMPHEP are interfaced to pythia for sub- 
sequent parton showering and hadronization. The MC 
events undergo a GEANT-based [15[ detector simulation 
and are reconstructed using the same algorithms as used 
for DO data. The effect of multiple pp interactions is 
included by overlaying data events from random beam 
crossings on simulated events. All MC samples except the 
W+jets are normalized to next-to-leading order (NLO) 
or next-to-NLO (NNLO) predictions for SM cross sec- 
tions; the tfjSingle t, and diboson cross sections are taken 
from Ref. [lj], Ref. 0, and the MCFM program [3, 
respectively. The Z+jets sample is normalized to the 
NNLO cross section [19(. The multijet background, in 
which a jet misidentified as an isolated lepton passes all 
selection requirements, is determined from data. In the 
muon channel, the multijet background is modeled with 
data events that fail the muon isolation requirements, 
but pass all other selections. In the electron channel, 
the multijet background is estimated using a data sam- 
ple containing events that pass loosened electron quality 
requirements, but fail the tight electron quality criteria. 
All multijet samples are corrected for contributions from 
processes modeled by MC. The multijet normalizations 
in the two lepton channels are determined from fits to the 
My distributions, in which the multijet and VU+jets rel- 
ative normalizations are allowed to float. The expected 
rate of multijet background is determined by this nor- 
malization, with an assigned uncertainty of 20%. 

Corrections are applied to the MC to account for dif- 
ferences from data in reconstruction and identification 
efficiencies of leptons and jets. Also, trigger efficiencies 
measured in data are applied to MC. The instantaneous 
luminosity profile and z position of the pp interaction 
vertex of each MC sample are adjusted to match those 
in data. The pr distribution of Z bosons is corrected 
at the generator level to reproduce dedicated measure- 
ments [20j. 

Other DO analyses of this final state apply additional 
corrections to improve the modeling of the VF+jets and 
Z+jets production in the MC Q. For the results pre- 
sented in this Letter, we choose not to apply those correc- 
tions in order to parallel the CDF analysis. We did, how- 
ever, study the effects of applying such corrections 21 1 
and find they do not alter our conclusions. 

We consider the effect of systematic uncertainties on 
both the normalization and the shape of dijet invari- 
ant mass distributions. Systematic effects are considered 
from a range of sources: the choice of renormalization 
and factorization scales, the ALPGEN parton-jet match- 
ing algorithm [22j, jet energy resolution, jet energy scale, 
and modeling of the underlying event and parton show- 
ering. Uncertainties on the choice of PDF, as well as 
uncertainties from object reconstruction and identifica- 
tion, are evaluated for all MC samples. 
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FIG. 1: (color online) Dijet invariant mass summed over elec- 
tron and muon channels after the fit without (a) and with (b) 
subtraction of SM contributions other than that from the SM 
diboson processes, along with the ±1 s.d. systematic uncer- 
tainty on all SM predictions. The % 2 fit probability, P(x 2 )> is 
based on the residuals using data and MC statistical uncer- 
tainties. Also shown is the relative size and shape for a model 
with a Gaussian resonance with a production cross section of 
4 pb at Mjj = 145 GeV/c 2 . 



In Fig. Q] we present the dijet invariant mass distri- 
bution after a fit of the sum of SM contributions to 
data. Other distributions are available in the supple- 
mentary material 21] . The fit minimizes a Poisson y 2 - 



function with respect to variations in the rates of individ- 
ual background sources and systematic uncertainties that 
may modify the predicted dijet invariant mass distribu- 
tion (23|. A Gaussian prior is used for each systematic 
uncertainty, including those on the normalization of each 
sample, but the cross sections for diboson and IF+jets 
production in the MC are floated with no constraint. The 
fit computes the optimal values of the systematic uncer- 
tainties, accounting for departures from the nominal pre- 
dictions by including a term in the fit function that sums 
the squared deviation of each systematic in units normal- 
ized by its ±1 s.d. Different uncertainties are assumed 
to be mutually independent, but those common to both 
lepton channels are treated as fully correlated. We per- 



TABLE I: Yields determined following a \ 2 fit to the data, 
as shown in Fig. [1] The total uncertainty includes the effect 
of correlations between the individual contributions as deter- 
mined using the covariance matrix. 





Electron channel 


Muon channel 


Dibosons 


434 ± 38 


304 ± 25 


W+jets 


5620 ± 500 


3850 ± 290 


Z+jets 


180 ± 42 


350 ± 60 


ti + single top 


600 ± 69 


363 ± 39 


Multijet 


932 ± 230 


151 ± 69 


Total predicted 


7770 ± 170 


5020 ± 130 


Data 


7763 


5026 



form fits to electron and muon selections simultaneously 
and then sum them to obtain the dijet invariant mass 
distributions shown in Fig. [TJ The measured yields after 
the fit are given in Table HI 

To probe for an excess similar to that observed by the 
CDF Collaboration [l|, we model a possible signal as a 
Gaussian resonance in the dijet invariant mass with an 
observed width corresponding to the expected resolution 
of the DO detector given by Ojj — aw- 
Here, aw^a and Mw 
the W 

11.7 GeV/c 2 and Mw->jj =81 GeV/c 2 from a simulation 
of WW — ¥ ivjj production. For a dijet invariant mass 
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jj resonance, determined to be aw-^jj 



resonance at Mjj = 145 GeV/c , the expected width is 
a-jj = 15.7 GeV/c 2 . 

We normalize the Gaussian model in the same way as 
reported in the CDF Letter We assume that any 
such excess comes from a particle X that decays to jets 
with 100% branching fraction. The acceptance for this 
hypothetical process (WX — > ivjj) is estimated from a 
MC simulation of W H — > Ivbb production. When testing 
the Gaussian signal with a mean of Mjj = 145 GeV/c 2 , 
the acceptance is taken from the WH — > ivbb simula- 
tion with Mh = 150 GeV/c 2 . This prescription is cho- 
sen to be consistent with the CDF analysis, which used 
a simulation of WH — > Ivbb production with Mh = 
150 GeV/c 2 to estimate the acceptance for the excess 
that they observes at Mjj = 144 GeV/c 2 . When probing 
other values of Mjj , we use the acceptance obtained for 
WH -> Ivbb MC events with M H = Mjj + 5 GeV/c 2 . 

We use this Gaussian model to derive upper limits 
on the cross section for a possible dijet resonance as a 
function of dijet invariant mass using the C'L S method 
with a negative log-likelihood ratio (LLR) test statis- 
tic [2^ | that is summed over all bins in the dijet invari- 
ant mass spectrum. Upper limits on cross section are 
calculated at the 95% confidence level (C.L.) for Gaus- 
sian signals with mean dijet invariant mass in the range 
110 < Mjj < 170 GeV/c 2 , in steps of 5 GeV/c 2 , allowing 
the cross sections for M^+jets production to float with no 
constraint. Other contributions are constrained by the 
a priori uncertainties on their rate, either derived from 
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FIG. 2: (color online) Upper limits on the cross section (in 
pb) at the 95% C.L. for a Gaussian signal in dijet invariant 
mass. Shown are the limit expected using the background 
prediction, the observed data, and the regions corresponding 
to a 1 s.d. and 2 s.d. fluctuation of the backgrounds. 



theory or subsidiary measurements. 

The Gaussian model is assigned systematic uncertain- 
ties affecting both the normalization and shape of the 
distribution derived from the systematic uncertainties 
on the diboson simulation. A fit |23[ of both the sig- 
nal+background and background-only hypotheses is per- 
formed for an ensemble of pseudo-experiments as well as 
for the data distribution. The results of the cross sec- 
tion upper limit calculation are shown in Fig. [5] and are 
summarized in Table [TT1 

In a further effort to evaluate the sensitivity for any 
excess of events of the type reported by the CDF Col- 
laboration, we perform a signal-injection test. We repeat 
the statistical analysis after injecting a Gaussian signal 
model, normalized to a cross section of 4 pb, into the 
DO data sample, thereby creating a mock "data" sample 
modeling the expected outcome with a signal present. 
The size and shape of the injected Gaussian model for 
Mjj = 145 GeV/c 2 relative to other data components is 
shown in Fig. [TJ 

The LLR metric provides a sensitive measure of model 
compatibility, providing information on both the rate and 
mass of any signal-like excess. We therefore study the 
LLR distributions obtained with actual data as well as 
the signal-injected mock data sample. The results of the 
LLR test in Fig. [3] show a striking difference between 
the two hypotheses, demonstrating that this analysis is 
sensitive to the purported excess. In the actual data, 
however, no significant evidence for an excess is observed. 

In Fig. |U we show as a function of cross section the 
p-value obtained by integrating the LLR distribution 
populated from pseudo-experiments drawn from the sig- 
nal+background hypothesis above the observed LLR, as- 
suming a Gaussian invariant mass distribution with a 
mean of Mjj = 145 GeV/c 2 . The p- value for a Gaussian 



FIG. 3: (color online) Log-likelihood ratio test statistic as a 
function of probed dijet mass. Shown are the expected LLR 
for the background prediction (dashed black) with regions 
corresponding to a 1 s.d. and 2 s.d. fluctuation of the back- 
grounds, for the signal+background prediction (dashed red), 
for the observed data (solid black), and for data with a dijet 
invariant mass resonance at 145 GeV/c 2 injected with a cross 
section of 4 pb (solid red). 
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FIG. 4: (color online) Distribution of p-values for the sig- 
nal+background hypothesis with a Gaussian signal with mean 
of Mjj = 145 GeV/c 2 as a function of hypothetical signal 
cross section (in pb). Shown are the p-values for the back- 
ground prediction (dashed black) with regions corresponding 
to a 1 s.d. and 2 s.d. fluctuation of the backgrounds and the 
observed data (solid black). 



signal with cross section of 4 pb is 8.0 x 10~ 6 , corre- 
sponding to a rejection of this signal cross section at a 
Gaussian equivalent of 4.3 s.d. We set a 95% C.L. upper 
limit of 1.9 pb on the production cross section of such a 
resonance. 

In summary, we have used 4.3 fb _1 of integrated lu- 
minosity collected with the DO detector to study the 
dijet invariant mass spectrum in events containing one 
W — > Iv (I = e or /it) boson decay and two high-pr 
jets. Utilizing a similar data selection as the CDF Col- 
laboration we find no evidence for anomalous, resonant 
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TABLE II: Expected and observed upper limits on the cross section (in pb) at the 95% C.L. for a dijet invariant mass resonance. 



M Si (GeV) 


110 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


170 


Expected: 


2.20 


2.01 


1.90 


1.78 


1.71 


1.64 


1.58 


1.52 


1.47 


1.40 


1.37 


1.31 


1.24 


Observed: 


2.57 


2.44 


2.35 


2.27 


2.19 


2.09 


2.00 


1.85 


1.69 


1.58 


1.46 


1.36 


1.28 



production of dijets in the mass range 110 — 170 GeV/c 2 . 
Using a simulation of WH — > Ivbb production to model 
acceptance and efficiency, we derive upper limits on the 
cross section for anomalous resonant dijet production. 
For Mjj = 145 GeV/c 2 , we set a 95% C.L. upper limit of 
1.9 pb on the cross section and we reject the hypothesis 
of a production cross section of 4 pb at the level of 4.3 s.d. 
In the case that the cross section reported by the CDF 
Collaboration is modified, we report in Fig. U the vari- 
ation of our p-value for exclusion of potential resonance 
cross sections other than 4 pb. 
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APPENDIX 

Fit of a 145 GeV/c 2 Dijet Resonance 

The DO data do not indicate the presence of a non-SM 
dijet resonance such as indicated by the CDF Collabora- 
tion. In the Letter we showed the fit of the SM predic- 
tions to the data. However, fitting only the SM contribu- 
tions could hide an excess if the systematic uncertainties 
allowed the SM contributions to be distorted in such a 
way that they filled in the excess. To study this question, 
we present a fit to the data of the SM predictions plus 
the Gaussian signal template with Mjj = 145 GeV/c 2 . 
The resulting dijet mass distribution is shown in Fig. [5] 
The fit is performed in the same way as described in 
the Letter (no constraint on diboson or VF+jets normal- 
izations), except that it now also includes the Gaussian 
signal template for Mjj = 145 GeV/c 2 with a freely float- 
ing normalization. The Gaussian model includes sys- 
tematic uncertainties affecting both the normalization 
and shape of the template, analogous to the systematic 
uncertainties for the diboson prediction. The best fit 
value for the Gaussian template yields a cross section of 
cr(pp — > WX) = 0.82^0 82 consistent with no excess. 
When we fix the diboson cross section to the SM predic- 
tion with a Gaussian prior of 7% on the rate, the best fit 
value for the Gaussian template yields a cross section of 
a(pp -> WX) = 0A2tolt pb. 

Kinematic Corrections to the Simulation 

The common tools used to simulate the predicted 
SM contributions perform well in general, but they 
have shortcomings. For example, different event gen- 
erators have different predictions for production angles 
and relative angles between jets in V^+jets and Z+jets 
events [Hj]. Thus, it is not unexpected that the simulated 
ly+jets and Z+jets samples do not perfectly model the 
angular distributions of jets. For analyses with looser 
selections, such as the search for WH production at 
DO these jet angular distributions show clear dis- 
crepancies between data and the simulated VF+jets and 
Z+jets events. Thus, these analyses use parameterized 
functions to correct the pseudorapidities of the two high- 
est pt jets and the A1Z — <J (A77) 2 + (A0) 2 separation 
between those two jets in W/Z+jets samples, and the 
transverse momentum of the W boson candidate, pr(W), 
in the VF+jets samples, to better model the data. 

The tight kinematic selection criteria employed in this 
analysis (e.g., pr{jj) > 40 GeV/c 2 ) remove much of 
the phase space in which the MC generators have diffi- 
culty modeling data (e.g., low pr{W)), greatly reducing 
the need for the kinematic corrections of the simulation. 
Therefore, the plots and results in the Letter do not use 
any of these kinematic corrections, which is consistent 
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FIG. 5: (color online) Dijet invariant mass summed over lep- 
ton channels when including the Gaussian model in the fit to 
data (a) without and (b) with subtraction of SM contribu- 
tions other than from the SM diboson processes, along with 
the ±1 s.d. systematic uncertainty on all SM predictions. 



with the CDF analysis. 

Although kinematic corrections are not required to 
achieve adequate modeling when applying the tight se- 
lection criteria of this analysis, modeling issues are prob- 
ably still present. In this section we present the results 
obtained when the kinematic corrections (derived from 
a selection similar to the search for WH production [3J) 
are applied to this analysis. 

The following figures are analogous to those in the Let- 
ter, except that the above mentioned kinematic correc- 
tions have been applied to the simulation. Figure |5] shows 
the dijet invariant mass distribution after the fit of the 
sum of SM predictions to data. The change relative to 
Fig. 1 in the Letter is not large, but improved modeling 
is evident in the higher x 2 probability. The resulting up- 
per limits on the cross section for production of a dijet 
invariant mass resonance are presented in Table IIIII and 
shown in Fig. [7] They are consistent with those in Fig. 2 
from the Letter. Figure [5] shows the LLR distributions 
analogous to Fig. 3 from the Letter. 

From this study we conclude that the kinematic cor- 
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FIG. 6: (color online) Dijet invariant mass summed over lepton channels after the fit (a) without and (b) with SM contributions 
subtraction other than from the SM diboson processes, along with the ±1 s.d. systematic uncertainty on all SM predictions. 
These distributions have the additional kinematic corrections applied to the MC. 



TABLE III: Expected and observed upper limits on the cross section (in pb) at the 95% C.L. for a dijet invariant mass resonance. 
These limits are derived with the additional kinematic corrections applied to the MC. 



M n (GeV) 


110 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


170 


Expected: 


2.35 


2.16 


2.05 


1.97 


1.88 


1.81 


1.73 


1.68 


1.65 


1.56 


1.52 


1.45 


1.42 


Observed: 


2.26 


2.02 


1.93 


1.83 


1.74 


1.64 


1.55 


1.48 


1.37 


1.27 


1.18 


1.09 


1.06 



rcctions to the simulation do not substantially change 
the result and we reach the same conclusion that there 
is no excess of dijet events in the DO data similar to that 
reported by the CDF collaboration. 

Additional Data-MC Comparisons 

Kinematic distributions presented in this Section are 
modeled without the additional corrections applied to the 



MC. Figure IH1 shows the dijet invariant mass distribution 
for the separate lepton channels after the simultaneous fit 
in these two distributions of the SM predictions to data. 
Figure 1 (a) in the Letter is the combination of these two 
plots. Figure [10] shows comparisons between data and 
simulation for other kinematic variables after the fit. 
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FIG. 7: (color online) Upper limits on the cross section (in 
pb) at the 95% C.L. for a Gaussian signal in dijet invari- 
ant mass. These results are derived with the additional 
kinematic corrections applied to the MC. Shown in the fig- 
ure are the limit expected using the background prediction, 
the observed data limit, and the regions corresponding to a 
1 s.d. and 2 s.d. fluctuation of the backgrounds. 
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FIG. 8: (color online) Log-likelihood ratio test statistic as 
a function of probed dijet invariant mass. These results are 
derived with the additional kinematic corrections applied to 
the MC. Shown are the expected LLR for the background 
prediction (dashed black) with regions corresponding to a 
1 s.d. and 2 s.d. fluctuation of the backgrounds, for the 
signal+background prediction (dashed red), for the data 
(solid black), and for data with a dijet mass resonance at 
145 GeV/c 2 injected using a cross section of 4 pb (solid red). 




FIG. 9: (color online) Dijet invariant mass distributions separately for the (a) electron and (b) muon channel after the 
simultaneous fit of these two distributions. 
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FIG. 10: (color online) Distributions of kinematic variables (combined electron and muon channels) evaluated using the results 
of a x 2 At of SM predictions to data for the dijet invariant mass distribution: (a) transverse W mass, (b) ATZ separation 
between the two selected jets, (c) lepton pr, (d) missing transverse energy, (e) highest jet pr, (f) second highest jet pr- The 
±1 s.d. systematic uncertainty on all SM predictions is presented by the cross-hatched area. 



